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Abstract 

We reformulate the matrix models of minimal superstrings as loop gas models on 
random surfaces. In the continuum limit, this leads to the identification of minimal 
superstrings with certain bosonic string theories, to all orders in the genus expansion. 
RR vertex operators arise as operators in a Z 2 twisted sector of the matter CFT. We 
show how the loop gas model implements the sum over spin structures expected from the 
continuum RNS formulation. Open string boundary conditions are also more transparent 
in this language. 
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1 Introduction 

String theories in low dimensions have recently come back into focus as simplified but still 
very instructive laboratories for string physics. The bosonic theories were originally solved 
through the double-scaling (=continuum) limit of hermitian matrix models, which implement 
a discretization of the bosonic worldsheet EQ-[S] The modern understanding of D-branes, 
open/closed duality, and tachyon condensation, together with crucial technical advances in 
Liouville CFT 0-^1], have inspired a reinterpretation of the double-scaled matrix model as 
the "open string field theory on an infinite number of decayed localized D-branes" |15| I16j. 
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In turn, this insight has led to the identification of the matrix models corresponding to non- 
critical superstring theories (c = 1 Type OB in [171 118j . c = 1 Type OA in 18 and c < 1 OA 
and OB in 19 j), and to a renewed interest in the subject (e.g. |2U|-|48|1. 

The Type OB cases are solved by unitary matrix models [SlEnil^lSlESlIS^ (equivalently 
hermitian models with a double-cut ansatz), while the Type OA cases are solved by complex 
matrix models [551 156| . This is surprising, as these matrix models are close cousins of the 
hermitian models that yield the bosonic string, and it is not at all clear in which sense (if any) 
unitary and complex models are providing a discretization of swperRiemann surfaces. This is 
not just a conceptual puzzle. These matrix models provide exact answers even in the presence 
of RR fluxes, while worldsheet CFT techniques would seem as problematic as in the usual 
critical string case. Connecting the matrix models with the continuum worldsheet formulation 
we can hope to learn something about the worldsheet realization of RR backgrounds. 

In this paper we make some progress in this direction. We focus on the (2, 4k) sequence 
of minimal superstrings, which are obtained by coupling (2, Ah) superminimal models to su- 
perLiouville theory and are related to one-matrix models [T15J - but we expect our approach 
to generalize to all minimal superstrings. The main technical novelty is to reformulate these 
matrix models in terms of a gas of self-avoiding fermionic loops on a random surface |57| . The 
random surface is the discretization of a bosonic worldsheet, with the loop gas as additional 
matter degrees of freedom. This implies that minimal superstrings are perturbatively equiv- 
alent to certain bosonic string theories, to each order in the genus expansion. The loop gas 
formulation allows to read off quite directly the matter CFT of these bosonic string theories. 
For example, in the simplest case of the (2, 4) model (pure supergravity, or c = 0) the matter 
theory is the free cm = —2 CFT of two Grassmann odd scalars Q + (z, z) and Q~(z, z). 

A general feature is that NSNS vertex operators map to vertex operators in the bosonic 
string with "integer" Liouville dressing; while RR vertex operators have "half-integer" Liou- 
ville dressings, their matter part belongs to a Z2 twisted sector of the matter CFT and they 
introduce cuts in the bosonic worldsheet, in analogy with the continuum RNS formulation. 
On higher genus surfaces, the loop gas model correctly implements the sum over spin struc- 
tures, again as expected from the RNS formulation. Thus, while we are formulating minimal 
superstrings (in the genus expansion) as certain bosonic string models, we seem to have got 
one step closer to the fermionic RNS formulation. It is an intriguing open question whether 
one can complete this program and recover our loop gas models from the point of view of 
discretizing surfaces with a metric and a gravitino. 

Open strings boundary conditions seem also more transparent in this framework. We 
find indications that the loop gas language accommodates naturally both "electric" and 
"magnetic" branes (77 = ±1 in the notations of |32|). while only the rj = — 1 branes are 
transparent in the one-matrix model |191 132j . We find that the 77 = — 1 branes are related to 
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Dirichlet boundary conditions for the fermionic matter degrees of freedom Q arising from 
the loop gas; they correspond to the standard resolvant of the matrix model. We propose 
that the rj = +1 branes are related instead to Neumann boundary conditions for the ® a 
fields, and correspond to a novel type of resolvant in the matrix model. 

For the closed string sector of the OA models, our claims reduce to well-established results 
|191 1581 QUI I45j . OA minimal superstrings with positive cosmological constant have no local 
RR vertex operators in the closed string sector ^]. The (2, 4k) OA theory is perturbatively 
the same as the minimal (2,2k — 1) bosonic string |191 1291 US*] . The OB models are obtained 
by taking a Z2 orbifold of the OA models (and vice- versa), where the Z2 acts as —1 (+1) on 
all RR (NSNS) vertex operators. Thus one may a priori expect the (2, 4k) OB model to be 
described perturbatively by a Z2 orbifold of the minimal (2, 2k — 1) bosonic string. The loop 
gas formulation makes this relation transparent, clarifying how exactly the orbifold works 
both in the matrix model and in the continuum worldsheet formulations. 

The detailed organization of the paper can be gleaned from the table of contents. Sections 
2, 3 and 4 deal with the closed string sector. Section 2 introduces the basic setup. Loop 
gas models for OA (OB) theories are deduced by gauge- fixing complex (two-cut hermitian) 
matrix models; the OB version is related to OA version by a Z2 orbifold that introduces 
a bicoloring of the random surface. Section 3 deals in detail with the c = case, where 
an interesting matrix model with Parisi-Sourlas supersymmetry EH EQ E3 E21 EH IM] 
plays an important role. Section 4 treats the continuum worldsheet formulation of minimal 
superstrings as bosonic strings. Section 5 focusses on D-branes, both in the discretized and 
in the continuum formalisms. Section El contains some concluding remarks. An appendix 
reviews properties of the c M = -2 "symplectic fermions" CFT ffifil IH7| 17)31 17HI I7T| 1721 173] . 
and details some of its uses in connection with minimal superstrings. 

2 The superstring from the loop gas 

We begin in this section by considering the Type OB (2, 4k) theories. We take the cosmological 
constant /x to be positive since our goal is to make contact with the continuum worldsheet 
formulation, and large positive fi is the usual perturbative regime. 1 These theories are solved 
by hermitian matrix models with a double-cut ansatz |19j . 

Recall that hermitian one-matrix model is defined by computing correlators with the 
measure 



1 The matrix model results suggest that large negative jj, corresponds to another perturbative regime, but 
it is not obvious how to reconstruct the corresponding continuum worldsheet theory. 




(2.1) 
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where M is an N X N hermitian matrix. The model is solvable by reducing the matrix 
integral to an iV-dimensional integral over the eigenvalues, which can then be evaluated 
using orthogonal polynomials theory (for a review see e.g. [73]). All correlators of basic 
U(N) invariant observables, e.g. TrM k and det(x — M), can be exactly evaluated this way. 
Alternatively, one can set up a perturbative expansion around a semiclassical vacuum. To 
expand around a one-cut vacuum (all eigenvalues in one potential well) it is convenient to 
shift M so that the bottom of the well is at M = 0, V'(M = 0) = 0. Treating cubic and 
higher terms of the potential as perturbations gives the usual expansion in ribbon Feynman 
diagrams, which can be organized according to their genus. Viewing these diagrams as 
discretized Riemann surfaces with a metric leads to the connection with minimal bosonic 
strings. 

Recently, the two-cut vacua of the hermitian one-matrix model for the left-right symmetric 
potential V(M) = V{—M) have been connected with Type OB minimal superstrings with 
fj, > |19j . The Z2 symmetry of the potential can be explained from the viewpoint of the 
open string field theory in the superstring. At present there is no connection between this 
matrix model and the viewpoint of discretizing superRiemann surfaces with a metric and a 
gravitino. 

2.1 Gauge fixing and the two-color loop gas 

To understand what kind of worldsheets this matrix model is discretizing, we now set-up a 
perturbative Feynman diagram expansion around the two-cut vacuum. It is natural to par- 
tially gauge-fix the action making M block-diagonal, M = diag(M^, Mr), and then expand 
independently in terms of Ml and M R around the two different critical points. 

The measure factor for the change of variables from M to Ml and Mr is standard (see 
e.g. |75| I76j.) The Jacobian determinant is the determinant of the linear action 

a^[a,dmg(M L ,M R )], (2.2) 

where a is a block off-diagonal matrix. Hence 

dM = det(M L (2) I Nr - I Nl ig> M R )det{M R ig> I Nl - I Nr ® M L ) dM L dM R , (2.3) 

where £ parametrizes u(N^)i!u(N R ) an< ^ ^ K ls ^ e K x K identity matrix. It is useful to 
rewrite this determinant as an integral over ghost variables, two Grassmann-odd Nl X N R 
matrices C\, C2 and their N R x Nl conjugate matrices Ci, 0%. The model ()2.1|) becomes 

J dM L dM R dCi dCi dC 2 dC 2 e~ S2 (2.4) 
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where 

S 2 = -Ti[V(M L ) + V{M R ) + CxMlCx - C 2 M L C 2 - M R C X C X + M R C 2 C 2 ] . (2.5) 
9 

Assume now that V(M) is symmetric under M — > — M and has a minimum at M = a and 
M = —a. It is convenient to define Ml = —a + <I>l and = a — $_r, 

M=f" a + ^ ° V (2.6) 

so that the two critical points correspond to 3>£ = and $k = 0. With this shift the ghosts 
acquire a kinetic term, 

S 2 = -Tr[V(<S> L -a) + V(<S> R -a)-2a(C 1 C 1 -C 2 C 2 ) (2.7) 
9 

+Ci$ L Ci - C 2 $ L C 2 + Qudd - <S> R C 2 C 2 ] . 

In analyzing this matrix model, we should remember that we are interested only in observables 
which are gauge-invariant with respect to the full U (TV) symmetry of the original one-matrix 
model for example 

TiM n = Tr(a - <5> R ) n + (-l) n Tr(a - $ L ) n . (2.8) 

If we add these operators to the action, we change the L and R potentials in a correlated way, 
corresponding to an analytic deformation of the original one-matrix model potential V{M). 
Keep also in mind that because of our choice of signs in ()2.6|) . operators which are even (odd) 
under $l <-> $r correspond to even (odd) powers of the original field M. 

The action (|2.7|) admits a good perturbative expansion. Feymann diagrams are easy 
to picture. Since the ghosts appear only quadratically, their propagators form self- avoiding 
closed loops in the Feymann diagram. There are no vertices mixing directly and <& R , 
which interact only through the coupling with the ghosts. The ghosts are Ni x N R matrices 
so in the usual double-line notation each ghost loop has two sides, on one side only <3?l legs 
are inserted, on the other only Hence any Feymann graph will consist of patches of two 
different "colors", L or R, separated by the ghost loops (see Figure^). The dual diagrams are 
random surfaces with the extra degrees of freedom of a gas of self-avoiding loops. There is a 
constraint on configurations, as there has to be an assignation of L or R color to each region 
bounded by a loop, with different colors on the two sides of each loop. In the evaluation 
of the diagram, each self-avoiding loop carries a weight of —2 (the minus sign because of 
odd-Grassmanality, and the factor of two because there are two copies i = 1,2 of the Cj, 



6 



M_L interaction vertex " M_L propagator 

M_R interaction vertex ' M_R propagator 

— — Ghost loop ^^^^^ Ghost propagator 

Figure 1: A Feymann diagram for the two-color 0{— 2) matrix model (|2,7|) . 



Ci) 2 . 

Loop gas models on random surfaces are well studied (e.g. jHZl El EE1 EO EH EH IHH 
1231 E3|)- I n the O(n) model one assigns a weight of n to each self-avoiding loop. Our case 
is then related to the 0(— 2) model. However the usual version of the 0(— 2) model |57] is 
formulated on a random surface with only "one-color". It is immediate to write the matrix 
model whose perturbative expansion generates the one-color 0(—2) loop gas [57] . 

J d^dCdC e~ Sl , (2.9) 

where 

Si = - Tr[F($ - a) - 2aCC + C<S>C + ^CC] (2.10) 

= i Tr[l/(A4) + CMC + MCC] . 
9 

(Here we have defined M = $ — a.) We obtain the two-color model (|2.7|) from (|2.10|) by 

2 It is clear that despite the sign differences in the action, the two flavors of ghosts play a symmetric role. A 
C2 ghost loop has a relative factor of (— ~Ly L + v R+ p with respect to a Ci ghost loop, where Vl and Vr are the 
numbers of L and R vertices in the ghost loop and P the number of ghost propagators. Since P = Vr + Vl, 
this factor is one. One could of course redefine variables to have Ci and C2 appear symmetrically in the 
action, but our choice is more natural in relation to the one-color model. 
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setting 

^(^ °), C-( ° Cl ). (2.11) 

The two-color model can be thought of as a Z2 orbifold of the one-color model. The orbifold 
operation introduces a twisted sector of operators odd under $l *-* At genus zero the 
untwisted sector (operators even under &l <-> Qr) is isomorphic to the one-color model, 
since "coloring" is trivial on the sphere. At higher genus the two theories differ also in the 
untwisted sector, since the two-color model has an additional constraint on configurations: 
the number of ghost loops linked to each non-trivial homology cycle of the surface must be 
even. 



2.2 Twist lines and RNS 

A familiar way to implement an orbifold is to relax the constraint and introduce "twist 
lines": On a surface of genus g, we pick 2g independent homology cycles Cj, j = 1, . . . , 2g, 
and write the partition function as an average over the 2 2g sectors labeled by the choices of 
signs Ej = ±1 for each homology cycle. Introducing into the partition function the factor 
or " Ylj e^ loops lmkmg c i implements the Z2 orbifold. 3 For example, for the genus one vacuum 
amplitude, we can write 

*»-. = £ (1 + ( " 1) T + ( " 1) ' > 42.ft - £ zfi(*. • (2.12) 

k,leZ fc,ZS2Z 

where the integers k and I denotes the number of ghost loops homotopic to the a- and 6-cycle of 
the torus, respectively; we also defined the amplitude Z^ 1 (fc, I) of the matrix model without 
any color constraint. The sector that includes the factor (— l) k , for example, corresponds to 
the torus with the twist-line along the 6-cycle. 

Matrix model operators consisting of even (odd) power of M are even (odd) with respect 
to the symmetry *-* (see (|2.8|) ) . In the formulation where we drop the color constraint 
and add twist lines, the appropriate signs for odd operators are accounted for by making 
them end-points for twist lines. Indeed a twist line between two insertion points will give a 
minus sign if they lie in patches of different color, a plus if they lie in patches of the same 
color (see Figure EJ). In the double-scaling limit of the two-cut matrix model, the even and 
odd operators correspond respectively to NSNS and RR vertex operators |19j . It is then 
compelling to identify the 2 2g sectors described above with the 2 2g choices of spin structures 
in the continuum RNS formulation. For example, the torus amplitude (|2.12|) can be regarded 

3 One says that the sector with ej — — 1 has a twist line along Cj, while the sector with ej = +1 does not 
have a twist line. 
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Figure 2: Twist lines and Ramond operators 



as that projected by the Type GSO projection 



Z g= i = Ttnsns 



1 + 

2 



1 + 

2 



(2.13) 



where fi,R is the left or right-moving world-sheet fermion number. Twist lines joining odd 
operators in the loop gas model clearly correspond to cuts joining spin fields on the RNS 
worldsheet. 

This seems a step forward in understanding how the matrix model is providing a discretiza- 
tion of a superRiemann surface. Clearly this program is not complete, since the two-color 
matrix model ()2.7|) is building bosonic random surfaces, with some extra matter degrees of 
freedom in the form of ghost loops. It would be very interesting to derive (|2.7|) from a super- 
geometric formulation that includes explicitly the path integral over the gravitinos. What we 
have learnt so far is that the minimal Type OB strings must have (at least perturbatively) a 
formulation as bosonic strings. 

2.3 Relation with OA 

With some hindsight, the perturbative equivalence of the OB theory to a bosonic string is not 
too surprising. The Z2 orbifold relating the two-color to the one-color model projects out all 
odd operators. In the continuum limit, this is the orbifold relating the OB to the OA minimal 
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superstrings, which projects out all RR closed string vertex operators. It is known that the 
string equations for OA (2, 4fc) model 4 are perturbatively equivalent to the string equations 
for the bosonic (2, 2k - 1) models (SHI El HH1 ES] • Hence we should expect that the (2,4k) 
OB theory is perturbatively equivalent to a Z2 orbifold of the bosonic (2, 2k — 1) model. 

This circle of ideas can be completed by showing directly that the complex matrix model 
for the OA theory is equivalent to the one-color loop gas. 5 The complex matrix model reads 

/ dMdM^-^ M ' M) , (2.14) 

where M. is a complex N x (N + q) matrix. Taking q = 0, we can gauge-fix to M. = 
hermitian and positive-definite. By a reasoning analogous to Q2.2l2.3j) . we are led to introduce 
the ghosts C, C and to reproduce precisely the action (|2.10|) . In deriving the one-color model 
by this route, the integration in (|2.10f) is restricted to positive-definite matrices M., 

f dMdCdC exp (--Tr[V(M) + CMC + MCC]), (2.15) 
JM>o \ 9 J 

where V(M) = V(M 2 ). This restriction makes no difference to each order in the perturbative 
expansion, and corresponds to the natural non-perturbative definition of the OA model. Let 
us mention for completeness the generalization to g 7^ 0. The gauge-fixing requires the 
introduction of some extra fermionic variables, the q x N matrix b and the N x q matrix b, 
leading to 

dM dCdC dbdb exp (-- (Tr[V(M) + CMC + MCC] + bM 2 b)) . (2.16) 



JM>0 \ 9 



The perturbative equivalence between the (2,4fe) OA with q = and the (2,2k — 1) 
bosonic models (5H1 El HH US] can be readily understood by establishing a map between the 
respective matrix models, as we now review. Diagonalizing M = diag(^j) and integrating 
out the ghosts in (|2.15|) . one has 

N fOC ! 

n j d »k no^ - ^? no^ + ^ e ~ n " k) ■ ^) 

k=l ^ i<j i,j 

The change of variables h{ = [if gives 

N 

TT / dh k T\(h i -h j ) 2 e-li V ^= [ dHe-- a V ^K (2.18) 
k=i Jo T<j Jh >° 

4 For q = 0, of course. More on the meaning of q momentarily. 
5 We thank N. Seiberg for explaining this to us. 
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This is almost the same as the hermitian one-matrix model for the matrix H = Ai 2 , with the 
differences that l)For generic potentials V(M), the action contains square roots of H; 2)The 
integration runs only over positive eigenvalues of H. However, for our purposes, we are only 
interested in even polynomial potentials V(M), which are polynomials in H. This is obvious 
from the point of view of the OA model, since V(y/~H) = V(H). Following the route from the 
OB theory, in going from the two-color to the one-color model we are restricting by definition 
to the "untwisted" sector of operators even under 3>£ «-> which map to even powers of 
M, TrAi 2n = TrH n . Moreover, the difference in the integration region has only effects which 
are non-perturbative in 1/N, and is thus immaterial to all orders in the genus expansion. 
Thus the one-color model is (perturbatively) equivalent to the standard one-matrix model. 
Its multicritical points 6 correspond in the continuum limit to the (2, 2k— 1) Virasoro minimal 
models coupled to Liouville, as anticipated. 

3 Loop gas models for pure supergravity 

We now elaborate on the simplest example. The (2, 4) OB theory (c = 0) is obtained from 
the two-cut matrix model with a quartic potential V(M) = \{M 2 — a 2 ) 2 [TJ2|. In the large 
N limit, for \i > one considers a two-cut ansatz for the eigenvalues. The critical point 
is reached when the two cuts meet at the origin with a zero of the resolvant R(z). In our 

4 

normalizations, this happens for gN — ► ^ from below. The double-scaling limit consists in 
setting gN = \ — efj,, e — ► 0, where /j, is the bulk cosmological constant, and zooming in 
around the origin of the z plane, z = ex. 

In the loop gas model (either in the two-color or one-color version), this critical point 
corresponds to the so-called "dilute phase". The loop degrees of freedom are becoming 
massless as we approach the critical point (the effective mass for the ghosts is the distance of 
the eigenvalue cut from the origin), so that the number of vertices occupied by the ghosts is 
diverging. Simultaneously the potential V is also becoming critical, which causes the number 
of unoccupied vertices to also diverge. (By contrast, the dense phase would correspond to 
the ghost degrees of freedom filling the random surface in the continuum limit.) In fact, this 
particular choice of quartic potential corresponds to a higher multicritical point, where both 
cubic and quartic vertices are diverging [^j. Standard formulas for the 0(n) model (see e.g. 
84 ) give the the central charge of the matter and Liouville continuum CFT, 

n = -2cos(^, Q = b+^, c L = 1 + 6Q 2 = 26 - c M , (3.1) 

6 multicritical polynomial potentials in H of degree k + 1, descending from multicritical even potentials 
V(M) of degree 2k + 2 in the original two-cut model 12. in 
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where b < 1 (6 > 1) in the dilute (dense) phase. 7 In our case we have n = —2 in the dilute 
phase, and the multicritical quartic potential corresponds to b = l/y2, cm = —2, cl = 28, 
which are indeed the values expected for the minimal (1,2) bosonic string theory. 8 Notice 
that we are landing on theory with b = l/y/2 rather than b = v2, so we are accurately 
referring to this theory as the (p, q) = (1,2) model (as opposed to (2,1)), in conventions 
where b 2 = p/q. This subtlety is immaterial for the closed string sector, which is our focus 
in this section, but plays a role for D-branes, as we shall see in section |SJ 

Following (|2.17I2.18|) . the one-color loop gas model with quartic potential V(M) can be 
mapped to a one-matrix integral with a gaussian potential, 

/ dHe-^ H - a ^. (3.2) 
Jh>o 

In this model, the large N eigenvalue density is the usual Wigner distribution, supported in 
the interval — (gN) 2 , \ + (gN)?]. As gN — > the left boundary of the cut reaches the 
origin (which is also the boundary of the integration region in (|3.2[) ). The double-scaling limit 
is again obtained by zooming in around the origin. The double-scaling limit of the gaussian 
model is well-known to yield the (2,1) minimal bosonic string [21|f)3], which is equivalent 
to topological gravity [S51 [HHl EH EE] , see HH for a recent simple derivation. 9 Our model 
is not quite the usual gaussian model in that we have a restriction of the integration region 
to positive H eigenvalues. We can think of having an infinite potential wall at H = that 
prevents eigenvalues to explore the region with H < 0. In the doubled-scaled variables, the 
wall is a finite distance away in units of \x and its effects are non-perturbative in g s . 

3.1 The supersymmetric matrix model 

The case of quartic potential V is especially interesting since we can view the one-color 
model (|2.1U[) as a matrix model with Parisi-Sourlas supersymmetry. Introducing the matrix 
superfield 

&(0,6) = $ + §C + 9C + 00F, (3.3) 
we can rewrite action Si in superspace [oTTl IfiT] . 

S 1 = - [ d6dOTr(d &d s $ + W($)), W(*) = -a* 2 + -* 3 . (3.4) 
g_J 3 

7 We are using conventions where the boundary cosmological constant operator is e 6 ^. The bulk cosmological 
constant operator is always e^-^ where b_ is the smaller of 6 and 1/6. Bulk operators must always obey the 
Seiberg bound, boundary operators need not. 

8 In principle n = —2 can be obtained by b 2 — l/(2m) for any integer m, but as shown in |79| the case of 
multicritical quartic potential corresponds to b = \/\/~2. 

9 The gaussian model is accurately labeled as (2,1) (as opposed to (1,2)). More in section 15.11 on the 
distinction between (1, 2) and (2, 1). 
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In fact, this super symmetric model is equivalent to the gaussian model (|3.2|) for any choice 
of superpotential W(<&), by the "Nicolai map" H = W'($) + a 2 (HI]. This formulat ion has 
also an immediate interpretation as a theory of matrix fields living in minus two dimensions 
|591l6Ull6T] . The perturbative expansion in superspace realizes the discretization of a Riemann 
surface embedded in minus two dimensions, which of course agrees with the value cm = —2 of 
the continuum string theory. The Feynman diagrams in superspace have a J dOjdQj integral 
at vertex j and a propagator 1/a 2 — (9j — 9k)(9j — Ok) /a. The propagators can be collected 
in an exponential like 

(3.5) 

This is the discretization of the continuum CFT of two free Grassmann odd scalars, with 
action 

— f cfze^d^d®^ (3.6) 



4-7T 

where a = ±. The superspace coordinates correspond to the zero modes of the 0° fields, 
6 = 9q, 9 = 9q . We refer to the appendix for more details on this CFT, which plays an 
important role in the following. 

We can now describe the continuum version of the Z2 orbifold that leads from the one- 
color to the two-color loop gas model. The effect of a twist line in the discretized theory is 
to flip the sign of the ghost propagators that cross it. This corresponds in the continuum to 
flipping the sign of the Q a across the twist line, that is, to performing the orbifold Q a — > —Q a . 
We conclude that (2, 4) minimal fermionic string is equivalent to a bosonic non-critical string 
defined by coupling to Liouville the matter CFT with cm = — 2 obtained by the Z2 orbifold 
of the Q a system. 



3.2 Open string field theory 

Finally we are in the position to point out a transparent open string field theory interpretation 
of the one- and two-color matrix models, in the same spirit as the proposal of |15l 1161 112 1 j 
for the c = 1 matrix model. Starting with the one-color model, we can regard ()3.4|) as the 
effective open string field theory on infinitely many localized D-branes of the (1,2) minimal 
bosonic string. We consider D-branes which have ZZ boundary conditions in the Liouville 
direction and Neumann boundary conditions for the Q a CFT. Since the coordinates 9, 9 that 
represents the positions of the D-brane are fermionic, we get the supermatrix model (|3.4[) . 
Incidentally, it is amusing to recover the matrix action as the action of Witten's cubic OSFT 
|89j (on ./V ZZ branes with Neumann boundary conditions for the Q a ), gauged-fixed to Siegel 
gauge and truncated to the lowest modes. Expanding the open string field as 

N 

I*) = E ( A '% c il% + C H °o c il% + Cij 0o"ci|O}ij + FijOpQ cilO),, + ...), (3.7) 
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one can evaluate the Witten action with the usual CFT rules 90. . The kinetic term cqLq^) 
gives precisely the quadratic terms in l|3.4|) . with the d$dg derivatives arising from the zero 
mode terms of Lq; while the interaction term ty, vl/) gives the cubic term of the superpo- 
tential. 

The Z2 projection (|2.11j) leading to the two-color model (|2.7j) can now be recognized as 
the usual rule jU^ to perform a Z2 orbifold on the D-brane worldvolume, 



where <r 3 = diag(I NL , -I Nr ). 

4 The continuum bosonic formulation 

The previous arguments show that we can view the (2, 4k) minimal superstrings as bosonic 
strings, to each order in perturbation theory. The OA theories map to the bosonic (2, 2k — 1) 
strings |19U58ll!Ml45| . while the OB theories to a Z2 orbifold of these models. We now analyze 
this correspondence in the continuum worldsheet formulation. 

4.1 The operator spectrum 

The first task is to establish a map between the physical vertex operators of the superstring 
and the physical vertex operators of the bosonic formulation. We shall see that the bosonic 
CFTs are not the usual {2,2k — 1) minimal models (at least not for the OB cases), but 
extensions with an infinite number of Virasoro primaries. A complete understanding will 
require the specification of the precise operator content and the construction of modular 
invariant partition functions. Here we consider the simpler problem of listing the Virasoro 
representations that appear in these CFTs. Our discussion partly overlaps with |19l I29j . In 
this subsection we assume k > 1 and come back to the special case k = 1 below. 

The primary operators of the (p, q) = (2, 4k) superminimal models are labeled in the 
usual Kac table notation as (r, s) = with 1 < s < 2k. Operators with odd s are in 

NSNS sector, while operators with even s are in the RR sector. The operator (1,2k) is the 
RR ground state, while (1, 1) is of course the identity. In the OA theory with fi > 0, only 
NSNS operators survive the GSO projection, while in the OB theory with \x > all operators 
are kept. To construct the physical tachyon vertex operators T rs one needs to gravitationally 
dress the superconformal primaries with a Liouville factor 10 e^ rs ^, with the exponent obeying 



10 Here we are using the same conventions as |32|. The background charge is given by Q = 6+1/6, b = \Jpjq 

and the Seiberg bound is at Q/2, for both the fermionic and models. The central charge of the Liouville sector 
is cl — 1 + 2Q 2 for superLiouville and cl — 1 + 6Q 2 for bosonic Liouville. 



*(0,0)=(T3-*(-0,-0)-<7 3 



(3.8) 
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the Seiberg bound (3 rs < Q/2. It is convenient to quote the rescaled dressings j rs = f3 rs /Q 
(the rescaling is such that the Seiberg bound is at 1/2). For the (2, 4fc) models, 

7l * = 272TTi)- 

In the bosonic (2,2k — 1) minimal models, operators are labeled as (f, 5) = (1,5) with 
1 < s < k — 1. Tachyon vertex operators T\g have rescaled Liouville dressings (again in units 
where the Seiberg bound is at 1/2) 

71S= 2^TT)- (42) 

We shall shortly need the more general formula 

2k + 1 - f(2k - 1) + 25 o , 

The matching between the bosonic and fermionic theories is now straightforward using KPZ 
scaling. KPZ scaling |92l I93j tells us that the partition function F g ({t n }) at genus g, as a 
function of the sources {t n } for the local operators of the theory, is a homogeneous function 
of degree Q(l — g), where each t n is assigned the weight 7 n equal to its rescaled Liouville 
(or superLiouville) dressing. Comparison of (|4.1|) and ()4.2j) leads to s = 25+ 1, i.e. to the 
identification 

T 12 ~ s +i = f ls . (4.4) 

As expected, the bosonic vertex operators map to NSNS vertex operators in the superstring. 
Notice that the operator Tn, which corresponds to the cosmological constant deformation of 
superLiouville theory, is absent in the minimal bosonic theory. It maps in the bosonic theory 
to (1,5 = 0), which is just outside the minimal Kac table. This is not surprising, since this 
operator is associated with a redundant deformation in the KdV hierarchy - which to each 
order in perturbation theory governs the OA theory as well. There are (k — 1) independent 
deformations in the integrable hierarchy, and k NSNS operators in the superstring; the A;-th 
operator (the one with lowest superLiouville dressing) must be redundant. In the bosonic 
theory this operator has a Liouville dressing equal to one half the Liouville dressing of the bulk 
cosmological constant (which has 5 = 1, of course) and it can be interpreted as a "boundary 
operator" [94 . In modern language, it corresponds to the boundary cosmological constant 
operator, which can be inserted on boundaries with FZZT boundary conditions. 

For the OA theory this is the end of the story. For the OB theory, we need to find the 
counterpart of the RR operators on the bosonic side. Inspection of (|4,1|) and (|4.3|) leads 
to match the superconformal labels (1, s) for s even with the bosonic labels (f = 2,5) with 
s = 2(5 — k + 1), i.e. to the identification 

T\2(~ s +i-k) — T25 , (4.5) 
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with k < s < 2k — 1. Clearly all of these operators are outside the minimal bosonic Kac 
table. Closure of the fusion rules will generically require to include all degenerate Virasoro 
representations (2, 5) and (1,5) with no restrictions on 5. In the OB case the enlargement of 
the Kac table is inevitable, and this suggests it may be natural to treat the OA case using 
a non-minimal table as well (see [55] for some relevant discussions). One can check that 
the bosonic fusion rules are compatible with the expected Z2 symmetry, with the operators 
(f , s) being even for f = 1 and odd for f = 2. It would be interesting to construct modular 
invariant partition functions for these non-minimal conformal field theories. 11 

4.2 c = from c = — 2 

The case of the bosonic (2, 1) model (more accurately (1, 2)), for which cm = —2, requires a 
separate treatment. This theory has been extensively studied (Ml Ell EH1 EM EDI IS E2 EH] as 
the simplest example of a logarithmic CFT. There is a free field representation of the model in 
terms of a pair of Grassmann odd scalars Q + (z, z) and _ (z, z). Modular invariant partition 
functions have been constructed, for various choices of operator content corresponding to 
various orbifolds of the Q a fields. Although these partition functions contain an infinite 
number of Virasoro primaries, they are rational with respect to a VF-algebra. It is conceivable 
that a similar story could be worked out for the other (2, 2k — 1) models. In the appendix 
we collect some useful formulas pertaining to the 0" system. 

The bosonic theory corresponding perturbatively to the (2, 4) OA model can be defined 
taking as matter CFT simply the free Q a system with periodic boundary conditions. The 
natural tachyon vertex operators are 

3 — n 1 _ 

O n = e *V n (dQ a ,d8 l3 )cc, n = 1,3,5,... (4.6) 

Here (p is the Liouville field and c, c the usual reparametrization ghosts of the bosonic string, 
and V n are matter primaries of dimension (n 2 — l)/8 built acting on the vacuum with Q 
oscillator. In the notation of the previous section, these primaries correspond to the (1,5) 
Virasoro representations, with s = (1 — n)/2. The primary 7 3 n= i is just the identity, its Kac 
label being (1,0) = (1,1) by the usual reflection property (f, 5) = (p — f, q — 5). 0\ is thus the 
bulk cosmological constant operator, and it maps to the only tachyon operator of the OA (2, 4) 
model, namely the cosmological constant operator of the superLiouville theory (s = 1 in the 
notations of the previous section). In this case the cosmological constant deformations of the 
fermionic and bosonic string theories are in correspondence with each other: the "conformal" 
background of the bosonic (1,2) model maps to "superconformal" background of the (2,4) 
model. This is not the case for the models with k > 1. 

11 There is some work on non- minimal c = theories, which are related to percolation |96l and the SU(2) 
WZW model at level zero |§T|. 
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For n > 1, the operators (|4.6|) are interpreted as "gravitational descendants" in the 
topological gravity language. From the viewpoint of the bosonic string theory, they are 
standard elements of the cohomology at ghost number (1,1). Gravitational descendants are 
often represented as elements of the cohomology at non-standard ghost numbers j^H], in fact 
this is necessarily the case if the matter CFT is truly minimal. In our case, the matter CFT 
has infinitely many primaries and we gain the ability to represent gravitational descendants 
as standard vertex operators of ghost number (1,1), which is perhaps a more convenient 
picture for actual calculations. This phenomenon is not new and has been discussed in the 
topological string theory literature in related contexts, see for example |99 l I100( HOlj . 

The theory corresponding to the (2, 4) OB model is obtained by the Z2 orbifold Q a — ► 
—@ a . The orbifold introduces a twisted sector built with Q a oscillators acting on a vacuum of 
dimension —1/8. We refer to the appendix for the expression of the torus partition function. 
The twisted sector can be decomposed in irreducible Virasoro representations of dimensions 
(n 2 — l)/8, where now n is even. All these Virasoro representations are also degenerate and 
correspond to the (0, s) elements of the Kac table, with s = —n/2. Tachyon vertex operators 
are obtained by dressing with the same Liouville exponential as in (|4.fij) . with n even. The 
tachyon with n = maps to the RR ground state of the (2, 4) model. Tachyons with n > 
are again interpreted as gravitational descendants. 



4.3 The torus path-integral 

As a check of our proposal, we can evaluate the torus partition function in the continuum 
bosonic formulation of the (2, 4) model and compare it with the matrix model result. The 
partition function in the continuum fermionic formulation of the (2, 4k) models has not yet 
been fully computed, because the odd spin structures are somewhat subtle. The results for 
the even spin structures are available |102j , and they match with the matrix model results JH] • 
For the OA theory with fi > 0, the spectrum has only NSNS operators, the torus partition 
function coincides with that of the bosonic (2, 2k — 1) theory, and we do not have anything 
new to add to the discussion in For the (2,4) OB theory with fj, > we should be able 
to compute the full torus partition function using our identification with the Z2 orbifold of 
the (1,2) model. The expected matrix model result is (for general k and [i > 0) 

2& = -^logt, (4.7) 

where t is the coupling of the operator of lowest dimension in the theory. 

This calculation is straightforward using the results of |102j . These authors computed the 
string theory torus path integral in the continuum formulation; in the case that the matter 



17 



CFT is a boson at radius R they found 



Z -Wv^)=-i(-! + ±f)lo g *. (4.8) 

As reviewed in the appendix, the Z2 orbifold of the Q a system has a partition function 
identical to a compact boson of radius R = v2c7- Thus the calculation of the torus path 
integral is almost identical to the one that leads to (|4.8|) with R = \f2ot' . The only difference 
is in the integral over the Liouville zero mode, which contributes the factor of — log t/(3, where 
(3 is the Liouville dressing of the most relevant matter operator. The torus path integral Z% 2 
of the Z2 orbifold of the (1,2) model can then be deduced from ()4.8j) if we account for the 
different values of (3. The most relevant operator of the (1,2) theory is the cosmological 
constant, which in our normalizations has (3 c =-2 = the most relevant operator for c = 1 
is again the cosmological constant, which has (3 c =i = 2. Hence 

Zz 2 = -^=L Z C =\V2) = -ilogt . (4.9) 

Pc=-2 8 

This agrees with the matrix model result (|4.T|) 12 for Z®^. 



5 D-branes 

A very interesting comparison involves the realization of open string boundary conditions in 
the loop gas model (and its continuum limit) versus the continuum fermionic formalism. The 
spectrum of FZZT branes in the continuum fermionic formulation OA and OB theories has 
been determined in [321 • Let us recall their results. In the (2,4) OB theory with [i > 0, there 
are three relevant FZZT boundary states |1()1 lll[ HI], labeled as (32] 13 

\x;t] = +1)b, \x;t] = -1,£ = ±1)b • (5.1) 

Here x is the boundary cosmological constant, which takes values in the auxiliary Riemann 
surfaces A42 4 introduced in JH2J. The label r\ = ±1 distinguishes the linear combination of 
left and right moving supercharges annihilating the boundary state ("electric" versus "mag- 
netic" branes) while £ = ±1 distinguishes a brane from its anti-brane. The \x;rj = +1) 

12 Curiously, the matrix model result 14.71 equals for all k (up to an overall factor of 1/4) the continuum calcu- 
lation of |1U2| for the minimal bosonic strings with modular invariant partition functions of type A p —iD 1+g /2, 
with the formal substitution p = 2, q = 2k — 1. This suggests a connection between the Z2 orbifold relating 
OA and OB and the Z2 orbifold relating the AA and AD bosonic theories, although the connection is not 
straightforward because here we are dealing with non-minimal CFTs. 

13 We are omitting the dependence on the label £ = sign(/i) of 32 : since we restrict our analysis to fi > 0, 
so that £ = 1, 77 = £77 = 77. 
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brane does not couple to the RR ground state and it is thus identical to its own anti-brane. 
Moreover, computations of one-point functions indicate that in BRST cohomology one can 
identify \x;r) = — 1; £) = | — x;rj = — 1; = £), and parametrize the independent rj = — 1 branes 
just by their value of x, forgetting the label £. Finally from the annulus amplitude, one finds 
that the open strings between branes with the same rj are in the NS sector, while open strings 
between branes with opposite rj are in the Ramond sector. Results for the higher (2, 4k) OB 
models with fi > are similar. A priori the boundary states would acquire labels corre- 
sponding to the matter primaries, however it is believed |32] (in analogy with the bosonic 
string) that in the BRST cohomology such labels are redundant and one can still parametrize 
the inequivalent boundary states byr/ = ±l,£ = ±l and x taking value in the appropriate 
Riemann surface -M-J^ik)' 

On the other hand, in the (2, 4) OA theory with \i > 0, the relevant boundary states are 

\x;r) = +l,£ = ±l) A , \x;rj = -l) A . (5.2) 

In the OA theory with \i > the RR ground state does not give rise to a local vertex operator 
[!§] . but it corresponds to the zero mode of the RR field, under which the rj = +1 brane is 
charged; changing the sign of £ reverses the charge of the brane. 



5.1 The rj = — 1 brane 

For both OA and OB, the FZZT branes with rj = —1 have been identified |19| I32j with the 
standard macroscopic loop operators of the corresponding matrix model. 



5.1.1 OA 

Let us consider the OA case first. The loop operator is 15 
WT~\xa) = "^Tr \og{M ] M-x A ) - -^Tr \og{M 2 -x A ) = -^Tr \og{H-x A ) , (5.3) 

where the arrow indicates the gauge- fixing procedure of section l2~3l In terms of the matrix H, 
this is the standard resolvant of the one-matrix model; perturbatively, results for the resolvant 
expectation values in (2, 4k) OA coincide with those in (2, 2k — 1) bosonic. However it would 
be wrong to attempt to match the auxiliary Riemann surfaces computed in [H^j between 
bosonic and fermionic models. This is because the results of 32] are valid for conformal 



14 However for the (2,4fc) OB models with k > 1, the rj = +1 brane does couple to RR fields and thus the 
distinction £ = ±1 is relevant. 

15 As usual, the operator creating a brane is actually the exponential of the macroscopic loop operator W(x). 
Also, by a slight abuse of notation we are using the same letter x in the matrix model operator W(x) and in the 
continuum boundary state \x; rj = —1), while the two are of course related by the standard renormalization. 
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Figure 3: Dirichlet boundary conditions for 0° in the loop gas model. The dashed lines 
represent the self-avoiding ghost loops. 

or superconformal backgrounds, but for general k the cosmological constant deformation of 
the bosonic theory does not correspond to the cosmological constant deformation of the OA 
theory. They do match for k = 1; indeed in this case the 77 = — 1, fj, > curve for OA coincides 
with the curve 2y 2 — 1 = x [22] of the (2, 1) bosonic theory (up to rescalings of x and y). 

It is clear that in our formulation in terms of the one-color loop gas, 1 (x) cuts 

a macroscopic hole which cannot intersect any of the self-avoiding ghost loops, since the 
matrices Ci do not appear in VF^ __1 (see Figure |3J). Let us now focus on the k = 1 case. In 
the supersymmetric formalism of section |3~T1 we can write 

Wr~\x A ) = ~Tr log((* - a) 2 - x A )\ e=§=0 . (5.4) 

In the continuum limit, this operator cuts a hole in the worldsheet with Dirichlet boundary 
conditions O a = for the cm = —2 CFT. Notice however that in the supersymmetric model, 
the more natural expression for the resolvant is 

W {1)2) (x) = ~Tr log((* -a- x)\ e=s=0 . (5.5) 

In the double-scaling limit, Wng) and 1 {xa) are equivalent with xa scaling as x 2 . This 

relation between x and xa reflects the fact that the continuum limit of the one-color gas 
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model gives what should be more accurately referred to as the (p, q) = (1, 2) minimal bosonic 
string, the one with b = \/p~Jq = l/v2; while the continuum limit of the gaussian model gives 
the (2, 1) theory with b = y/2. There is no distinction between (1, 2) and (2, 1) in the closed 
string sector but open string boundary conditions are labeled differently. The boundary 
cosmological fis constant of the (1,2) model is equal to the dual boundary cosmological fis 
of the (2, 1) model, and vice-versa. Now Mb'' = Mb = 2(mb' 2 ^) 2 ~ M> which is compatible 

(2 1) (1 2) 

with the identifications xa ~ Mb an d x ~ Mb • FZZT branes with Dirichlet boundary 
conditions for Q a in the (1,2) have been considered in |34| . where is is shown that the full 
cubic OSFT on them localizes to the Kontsevich model. The authors of |44j reproduced this 
result starting with the resolvant ()5.3|) for the gaussian matrix model in terms of the matrix 
H, and performing the double-scaling limit; they observed that in their procedure one lands 
to the (2, 1) model with b = \/2, a fact that we have just explained from the relation between 
the gaussian and loop-gas matrix models. 16 

More generally the full superfield loop operator, 

W (h2) (x,p a ) =-±Trlog((&-a-x)\ e=Pi§= - p , (5.6) 

gives in the continuum limit the FZZT brane tensored with the Dirichlet boundary state 
|Dir(p a )) for the Q a system (see (|A.13|) in the appendix). 

5.1.2 OB 

Let us now consider the 77 = — 1 brane for the OB model. This is naturally identified with the 
macroscopic loop operator of the two-cut matrix model, 

WT _1 (x B ) = -^Tr\og(M -x B ) = --Trlog(-XB-a + ^ L )--Tr\og(-XB + a-^ R ) 

N ^l{a + XB) N ^ I (a — xb) N 

Again in the two-color loop gas model the ghost loops do not intersect the boundary (Figure 
I3J). It is also clear from Q5.7|) that the expansion of W b ^~ 1 (xb) contains both even and odd 
operators under <J>^ «-> which in the continuum limit are identified as NSNS and RR 
operators. This is just as expected for the rj = — 1 branes, which couple both to NSNS and 
RR closed operators (see (|5.1j) ). We can be more precise. In the expansion of W b ^~~ 1 {xb)-, 

16 The correspondence between the gaussian matrix model and the (2, 1) string theory can also be deduced 
using the general loop gas formula 13.11 . Besides the dilute phase realization with n = —2 and b — l/v2, there 
is a dense phase realization of cm = —2, with n = and b — Since n = 0, we have just the one-matrix 
model, and since we are in the dense phase, the potential is not at a critical point - and we may as well take 
it to be gaussian without changing universality class. 
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we weight a boundary of length I differently according to its color, either with a factor 
proportional to l/(a + xf for the color L or with a factor proportional to l/(a — x) 1 for the 
color R. To express this fact in the formalism of the unconstrained loop gas with twist lines, 
we must extend to sum to new sectors, allowing for twist lines that end on a boundary. The 
weight for a boundary of length I must be taken 



1{xb + o) 1 I(xb-cl) 1 

where e = +1 in the sector where no twist line ends on the boundary and e = —1 if it does. 
This fits perfectly with the continuum RNS formulation, where the rj = — 1 branes have both 
an RR and an NSNS term in their boundary state and thus one must sum over the different 
sectors in the closed channel. If we remove the distinction between L and R (that is, project 
out the odd operators), 

WT~\xb) — -^Tr log((* - a) 2 - x\) = WT~ l {x A = 4) • (5-9) 

The relation 

X A — IS clS expected pil32]. 

We can now analyze the OB rj = — 1 brane in the continuum bosonic formulation. The 
boundary state is obtained by taking the Z2 orbifold of the FZZT brane with Dirichlet bound- 
ary conditions Q a = 0. The orbifold operation yields two fractional branes with opposite 
coupling to the twisted sector, which we identify with the \rj = — 1,£ = ±1)b branes. This 
clearly agrees with the discretized expression (|5.7|) . which cuts holes in the worldsheet with 
Dirichlet b.c. for the Q a , and on which twist lines can end. 



5.2 The rj = +1 brane 

The brane with rj = +1 has so far been more puzzling, since it does not appear to correspond 
to standard macroscopic loop operators in the one- or two-matrix models |191 132j . In our 
formalism (focussing on the (2, 4) case) there is a natural guess for what this brane should 
correspond to, namely the other natural choice of boundary conditions in the Q a system: 
Neumann boundary conditions. The most compelling piece of evidence for this proposal 
comes from consideration of the open string spectrum between Neumann and Dirichlet branes 
for the Q a system, which lies in the twisted sector (see (jA.16|) ) - and hence has "half- integer" 
Liouville dressing (i.e. n even in the notations of section 4.2). This agrees with the fact that 
the open strings between rj = — 1 and rj = +1 branes are in the Ramond sector. Since the 
identification of the 77 = —1 brane with Dirichlet b.c. for Q a is beyond any doubt, we seem 
to be forced to identify the r\ = +1 brane with Neumann b.c. for Q a in order to reproduce 
the correct open string spectrum. 
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In the matrix model, we can find an expression for the Neumann macroscopic loop op- 
erator by appealing to the string theory intuition. We regard the matrix model as the open 
string field theory on Neumann ZZ branes as in |15l I16j . To a Neumann FZZT brane we 
add the open strings stretched between the Neumann FZZT brane and the N Neumann ZZ 
branes. Such open strings are expected to be fermionic |44| . and since we have Q Neumann 
b.c. on both branes, they will depend on the Q a zero modes. Thus we need to introduce a 
Grassmann odd open string field ir(8, 9) in the fundamental of U{N). In the supersymmetric 
formulation we think of it as a new superfield. We are led to the matrix model 

J d$dir e ~( s *+ s *) (5.10) 

with 

= J dOS Tr(d g $d s G> - a* 2 + * 3 ) , = J d6d6 (d e ^d 5 7r + wnK + tt^tt) , (5.11) 

where w is a mass parameter. In the Feynman diagram expansion, the new field introduces 
boundaries which are free to fluctuate in the Q a directions. In the loop gas picture, now 
there will be ghost lines attaching to the boundary. The tt superfield can be integrated out, 

J d$ • det($ + w)~ 2 ■ det f ($ + w) 2 - C-^^C - F^j ■ e~ 5 *. (5.12) 

One can imagine to further integrate out F, C and C. This would lead to a complicated 
expression involving a sum of multitrace combinations of the lowest-component matrix <£. We 
leave a detailed study of this new resolvant and of its double-scaling limit for future work. 

It would be nice to understand in the framework of the continuum bosonic formulation 
why Neumann brane for Q a should be a source of RR charge in the OA model, but not in the 
OB model (see (j5.lU5.2jl ). At present we have no crisp answer. RR charge is related to the very 
subtle RR ground state, which both in the supersymmetric and in the bosonic formulations 
precisely saturates the Seiberg bound. A logical possibility is that while minimal superstrings 
are perturbatively equivalent to bosonic strings in the closed string sector, this equivalence 
breaks down for D-branes, which are non-perturbative objects. This is a valid objection 
for ZZ branes. 17 In the matrix model ZZ branes are related to stationary points of the 
effective eigenvalue potential, and the non-perturbative restriction on the range of eigenvalue 
integration will affect their spectrum. However, we expect that FZZT branes should have a 
good description in the bosonic formulation. Intuitively this is because FZZT branes can be 
continuously connected to the closed string vacuum by sending the boundary cosmological 
constant to infinity, and in that limit they admit an expansion in closed string operators. We 
hope that a better understanding of the matrix model (|5.11|) and of its continuum limit will 
clarify this issue. 

17 We thank N. Seiberg for raising this point. 
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5.3 Comparison with previous work 

We would like to make contact with previous discussions of boundary conditions for loop gas 
models (E2l ESI EI] • Loop gas models can be formulated as "height" models (see [EU for a 
recent discussion), where configurations are specified by assigning a height variable to each 
vertex of a triangulated random surface (the heights of two neighboring vertices differing at 
most by one unit). The self-avoiding loops are the domain walls separating regions of the 
surface of different heights. The height model of [Hi\ is rich enough to describe the (dicretized 
version of) non-minimal matter of central charge cm < 1 plus Liouville CFT of central charge 
cl = 26 — cm- In the continuum limit, the height variable of 83] goes over an anomalous 
free boson X which provides a Coulomb gas representation of the matter CFT. In the height 
model, the natural boundary conditions are "Dirichlet" (fixed height at the boundary) and 
"Neumann" (free sum over different heights at the boundary). "Dirichlet" boundary con- 
ditions are equivalent to not allowing the self-avoiding loops to touch the boundary. In the 
continuum limit, they correspond to Dirichlet boundary conditions for the matter field X, and 
FZZT boundary conditions for the Liouville BCFT. The "Neumann" boundary conditions of 
|82 [ 1831 IS3] correspond instead to summing over all possible configurations, with self-avoiding 
loops allowed to end on the boundary. In the continuum limit, these boundary conditions go 
over Neumann boundary conditions for the Coulomb gas field X, and dual FZZT boundary 
conditions for the Liouville [S3]. Dual FZZT boundary conditions mean that the dual defor- 
mation e^l h is turned on the worldsheet boundary as opposed to e b ^. 18 In fact, we know from 
the FZZT boundary CFT 10 that both the e b ^ and e^l h deformations are actually turned 
on in the exact solution, so in the Liouville sector the distinction between ordinary and dual 
FZZT b.c. is immaterial (up to relabeling parameters); however in the matter sector, the 
distinction between "Neumann" and "Dirichlet" b.c. for X is significant. 

The boundary conditions that we have identified with the rj = — 1 FZZT brane, namely 
Dirichlet boundary conditions for Q a , are clearly the same as the "Dirichlet" boundary con- 
ditions of [S3]. However it is not at all obvious that Neumann boundary conditions for the 
0° system (and their discretized version l|5.11|) ) are the same as the "Neumann" bound- 
ary conditions of 84 . The continuum limit of the height model of [H3J (specializing their 
formulas to n = —2 in the dilute phase) gives a free anomalous boson X with cm = —2 
coupled to Liouville. Clearly "Neumann" boundary conditions for X have nothing to do 
with Neumann boundary conditions for Q a . However this comparison with the formulas of 
|84j is naive, because their model is different from ours in some crucial ways, even in the 
one-color case. The super symmetric matrix model corresponds to a 0{— 2) loop gas model 

18 Recall that our conventions are such that b 2 = p/q (so b can be either smaller or bigger than one according 
to the model), the boundary cosmological constant operator is always e 6< * and the dual boundary cosmological 
constant operator e^^ b . In |84| different conventions are used, with b < 1 by definition. 
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with multicritical quartic potential |79| ; whereas the model of [S3] is based on a triangulation 
of the surface. Another (probably related) point is that the Q a CFT does not need screening 
charges, whereas the X CFT does. It will be very interesting to extend the formalism of |84| 
to our case. It is conceivable that a suitable modification of the loop equation approach of 
[821 184j to "Neumann" boundary conditions may turn out to be equivalent to Q5.11JI . and 
provide a way to do concrete calculations. 

6 Outlook 

We have seen that the matrix models of minimal superstrings have a very natural formulation 
as loop gas models on random surfaces. This language has many virtues. 

First, it clarifies the Z2 orbifold relation between the OA and OB matrix models, which 
amounts to going from an ordinary random surface to a bicolored random surface. Second, it 
gives in the continuum limit a formulation of minimal superstrings as bosonic string theories, 
to each order in the perturbative expansion. The main news here are for the OB cases, since 
the OA cases have no local RR operators in the closed string spectrum and have been known to 
be perturbatively equivalent to minimal bosonic theories j^l ESI HH1 133 • I n the OB examples, 
RR operators are in the twisted sector of the Z2 orbifold and we have seen very explicitly how 
the loop gas model implements the sum over spin structures; in the simplest case of c = 
this is the sum over spin structures of the Q a system. One of the most intriguing directions 
for future work is to see if one can recover the loop gas model from a direct discretization 
of the RNS worldsheet. The @ a fields seem vaguely reminiscent of gravitinos, in that RR 
operators make them multi-valued - but they have, of course, the wrong spin. 

Finally, and this may turn out to be the most practical aspect of our work, loop gas 
models appear to give a more transparent description of open string boundary conditions. 
We find it very plausible that the r/ = +1 branes, which are mysterious in other formulations, 
correspond to Neumann boundary conditions for the Q a system. Further work is needed to 
completely settle this issue. A promising direction is a detailed study the new macroscopic 
loop operator that we have proposed in section HT21 

Several extensions of the results of this paper can be contemplated. More general (p, q) 
minimal superstrings should be described by two-matrix models, and their gauge-fixing should 
lead to loop gas models. It would also be nice to see if a similar gauge-fixing procedure 
as described in this paper can give some insight into c = 1 string theory. Some progress 
in connecting the c = 1 Type OA theory with the RNS formulation has been made in the 
interesting paper [3S], and it would be nice to understand the relation between their approach 
and ours. 

The theme that underlies the whole subject is open/closed duality. The doubled scaled 
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matrix models dual to non-critical strings are interpreted as the effective open string field 
theories on a large number of decayed ZZ branes. In this paper we have seen a new sharp 
application of this idea. The one- and two-color loop gas models are recognized as the 
effective worldvolume theories on the appropriate ZZ branes; it is satisfying to see that the 
Z2 orbifold relating them is the usual string theory jHJ. The other class of Liouville branes, 
namely the FZZT branes, offer an alternative route to open/closed duality [HI]) which we 
have not explored in this paper. The OSFT on the FZZT branes of (1,2) minimal bosonic 
string localizes to the Kontsevich matrix model |34| I44j . The speculation of (Hlj is that 
"topological" matrix models a la Kontsevich are generically related to FZZT branes. In the 
context of minimal superstrings this raises many natural questions. We expect that the study 
of FZZT branes in minimal superstrings will lead to interesting generalizations of topological 
matrix models. 

As usual with solvable models, the most important question is which features of the exact 
solution generalize to more physical situations. Open/closed duality is certainly one such 
feature. Other lessons specific to minimal superstrings should be sought in the study of 
exact RR backgrounds. In the critical RNS string, there are two related difficulties with RR 
backgrounds: they introduce cuts on the worldsheets and thus cannot be exponentiated in 
any obvious fashion; they break superconformal symmetry and thus the very rules for string 
theory computations are unclear. In this paper we have seen that for minimal superstrings 
the second difficulty can be circumvented by going to an equivalent bosonic formulation. This 
simplification may well be an artifact of the simplicity of these models. However, the first 
problem is still very much with us in the continuum bosonic formulation. Remarkably, the 
matrix model manages to produce an exact answer. A paradigmatic example is the Ising 
model with external magnetic field. Either as a continuum CFT, or as a spin system on 
a regular lattice, this is a notoriously difficult problem. The model becomes vastly simpler 
once formulated on a random lattice 103j; it can be easily mapped to a two-matrix model 
with asymmetric potential, and solved exactly. Summing over triangulations before summing 
over the spin degrees of freedom is the winning strategy. We can only speculate about the 
implications for critical superstrings. Perhaps the message here is that in trying to formulate 
RR backgrounds on a fixed Riemann surface, we are addressing a more difficult problem than 
we really need to solve. 
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A Appendix: The @ a CFT 

In this appendix we review some basic facts about the "sympletic fermions" CFT |661 1671 1681 
1691 17U1 1711 1721 173j and elaborate on its connections with minimal superstrings. 

A.l Bulk theory 

Q + (z, z) and Q~(z,z) are non-chiral, Grassmann odd fields of dimension zero, governed by 
the free action 

— j d 2 ze aP d@ a d@ p . (A.l) 

Here a, (3 = ±, e ±=F = ±1, e^e^ 7 = <5J. Hermitian conjugation acts as (0 + )^ = 0~. This 
CFT differs from the more familiar r]£ system in the treatment of the zero modes. One 
can identify £(z) = @ + (z,z), rj(z) = dQ~(z,z) and 77(f) = d@~(z,z), but there is a 

mismatch in the zero mode sector. While the r]£ system has two chiral zero modes (one for 
the chiral field £ and one for the antichiral field £), the Q a system has two non-chiral zero 
modes, one for + and one for 0~. 

Let us first consider the "untwisted" theory, with periodic fields Q a (z, z) = @ a (e 2m z, e~ 2m z). 
The mode expansion reads 

OO / CX. — Q \ 

Q a (z,z) = 9% + X " ln\zf + i £ fe + ^j. (A.2) 

n=-oon^0 

The modes obey the obvious canonical anticommutation relations 

{xS,eg} = <?P, {x^,X^} = me^S m+ni o, xi} = ™e^5 n+mfi . (A.3) 

Notice that in correspondence with the two non-chiral zero modes 8q, there are two canon- 
ically conjugate non-chiral momenta Xoi *- e - one has the identification X q = Xo, which is 
required by locality. The Fock space is obtained by acting with the creation operators X - n 
and X-n ( n > 0); on the space of ground states, which is spanned by 

in), \e a ) = e^\n), | w ) = ^o|n>. (a.4) 
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Here \Q) denotes the SL(2,C) vacuum. This system provides the simplest realization of a 
logarithmic CFT. The states |f2) and \uj) form a logaritmic pair, spanning a two dimensional 
Jordan block for Lq, 



Lq\u) =Lq\u) = Lo\Sl) = L„|fi) = 0. 



(A.5) 



Let us also record the inner products 



(n\n) = (9 a \n) = (e a \uj) = o, (n\u) = i, {e a \d p ) = e aP , (u\u) = k. (a.6) 



Here k is an arbitrary normalization constant, corresponding to the freedom to redefine 
\lo) — ► \lo) + a|fi). If we stick to the definition in (|A.4|) . then k = 0. 

The path-integral is zero unless one saturates the fermionic zero modes. The simplest 
non-zero correlators involve a single insertion of the logarithmic identity w, 



where J-k(wk,Wk) is an arbitrary local operator containing no zero modes. Such a correlator 
does not depend on the position of the ui insertion, since taking a derivative with respect 
to z gives a correlator that has unbalanced zero modes and is thus zero. On the sphere, 
a correlator with a single uj insertion factors into a holomorphic rational function of the 
coordinates Wi times a antiholomorphic rational function of Wi. Correlators involving more 
than one u insertion contain logarithmic terms log \ wi — Wj\, but they will not be important 
for us. Notice that even on higher genus Riemann surfaces one insertion of u is sufficient to 
saturate the zero modes. Moreover, on a surface of arbitrary genus, the path integral in the 
non-zero mode sector will vanish unless the number of a = + insertions equals the number 
of q = — insertions. This is in contrast with the r/£ system, where on a surface on genus g 
one needs to saturate one chiral £ zero mode (the constant mode) and g chiral ij zero modes 
(corresponding to the g holomorphic one- forms); and similarly one £ and g f) antichiral zero 
modes. The torus partition function with periodic boundary conditions and one insertion of 
the logarithmic identity reads 



which is just the inverse of the partition function for two free bosons, and is clearly modular 
invariant. 

For our purposes, the Q a theory with periodic fields is the matter CFT with cm = —2 
that enters in the construction of the (1,2) bosonic string. As usual, the string theory is 
defined by coupling the matter CFT with Liouville CFT of cl = 26 — cm = 28 and the be 




(A.7) 
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ghosts of cq = —26. As observed long ago in )87) . the connection with topological gravity is 
provided by taking the Liouville CFT to have /x = and performing the bosonization 19 

/3 = «9e + e^, 7 = dG~e~^, (A.9) 

where (ft is the Liouville field. The bosonic ghosts /?7 have central charge 26 and are the 
superpartners of the fermionic ghosts be in the topological gravity multiplet of 85 . The 
bosonization (|A.9j) does not involve the @ a zero modes. If we want the path integral over 
<p and Q a to reproduce the /?7 path integral, we should simply insert a single logarithmic 
identity u in every correlator, as in ()A.7|) . Similarly, since the Liouville path integral for fi = 
is proportional to the (infinite) volume Vol(^o) of the Liouville zero mode, in making contact 
with the /?7 system we need to divide out by Vol(0o). 20 With these rules, the simplest 
non-vanishing correlator is the sphere amplitude of three cosmological constant operators 
0\ = cce^ 2 ^ (with the implicit extra insertion of u> at some arbitrary point on the surface). 
This correlator is of course just a constant, in accord with the fact that the topological gravity 
|86| partition function at genus zero is proportional /i 3 . 

The Z2 orbifold Q a —* — Q a projects out all odd states and introduces the twisted sector 
with antiperiodic boundary conditions @ a (z,z) = —Q a (e 2m z,e~ 2m z). There is a unique 
twisted sector ground state \RR), with dimension (—1/8,-1/8). The full torus partition 
function of the Z2 orbifold, obtained by summing over spin structures, can be written as |66| 



,00 _ _ _ _ 

-J— y q l(V2m + n/V2f-±(V2m-n/V2)^ (AJQ) 

7?(r)r ^ 

1 / v m,n=— 00 

which one recognizes as identical to the partition function of a compact boson at radius 
R = V2a'. The Z2 orbifold theory is identical to the "triplet model" of |69[ I7D| , which was 
originally introduced by an algebraic construction involving a triplet of W3 generators. In 
the explicit free field representation provided by the symplectic fermions, the W3 generators 
are the spin 3 operators <9 2 e+<9G+, (d 2 Q + dO~ + d 2 Q-d& + ) and d 2 Q-d®-. The OPE of 
these operators close with the stress tensor to form an enlarged chiral algebra. The theory is 
rational with respect to this extended algebra. 

A. 2 Boundary states 

Open string boundary conditions for symplectic fermions have been studied in [7T1 1721 173] . 
Boundary states preserving the VF-algebra are either Neumann or Dirichlet, 

(x^-X-J|Dir)=0, { X a m + x_J|Neu)=0, Vm > , (A.ll) 



19 The arguments of |87| are actually phrased in terms of 77^ system. 

20 The discrete analog of dividing out by Vol(0o) is to divide out by the logarithmic divergence | log A| found 
in the matrix model (where A is the bare cosmological constant), as is done in 1551 . 
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with the understanding that Xo = Xo • I n the untwisted theory, these conditions are solved 
by 

|Dir) = N D ^ V [Y;\xZ k xt k + \xZ k xt k )\U) D , (A.12) 
\fc>o " / 

|Neu) = Mn exp ^ ~\xZ k xZ k - rXZ k xZ k . 

\fe>o ' / 

Here |0)d indicates any of the vacua in (|A.4D . since for the Dirichlet case the m = condition 
in l|A.ll|) is trivially implied by locality (xo — Xo)- ^ n ^ ac ^ it is natural to assemble the four 
possible Dirichlet boundary states into the linear combination 

|Dir(p a )) = N D exp(e a pp a X P) exp I ]T ~xZ k X± k + \xZ k xt k I M , (A.13) 

\A;>0 ' / 

where p Q are Grassmann numbers. Since #g|Dir(/j a )) = p a |Dir(p a )), this is interpreted as 
the Dirichlet brane at "position" Q a = p a . In the Z2 orbifold theory, the Dirichlet branes 
built on \9 a ) are projected out. Two new states satisfying (jA. 1 If) are built on the twisted 
sector ground state, 

ex P (j2^xZ k xt k ± \xZ k xz}j \RR) , (A.14) 

where the plus (minus) signs correspond to Dirichlet (Neumann) boundary conditions. (Here 
X k , X k are half- integer moded.) The consistent boundary states obeying Cardy conditions 
are specific linear combinations of the states in (|A.12|) and in (|A.14|) . we refer to |72l I73| 
for details. (The construction of the boundary states is somewhat more subtle than in the 
paradigmatic minimal model example, because of the non-trivial Jordan block structure of 
some of the representations, but can nevertheless be carried out). 

The essential point to emphasize here is that the open string spectrum between a D-brane 
with Dirichlet boundary conditions and a D-brane with Neumann boundary conditions lies 
in the twisted sector. Consider the cylinder amplitude 

(Neu\(qVY 0+Lo+1/6 \VHp a )) = ^ TT( 1 + f 1 ) 3 = ^ • (A.15) 

4 ± \ 2 ri{T) 

Modular transformation to the open channel gives (with NdNn = 2) 

MlI = q c/^ q -v8 fj (1 _ ^+1/2)2 . (A . 16) 

n=0 

Here c = —2. We recognize the spectrum of open strings obtained acting with half-moded 
oscillators {x°L k \ a = ±, k = 1/2, 3/2, . . .} on the twisted vacuum. 



30 



References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 

[9 

[io; 

[12 
[13 
[14 

[15 
[16 

[17 
[18 



M. R. Douglas and S. H. Shenker, Nucl. Phys. B 335, 635 (1990). 

D. J. Gross and A. A. Migdal, Phys. Rev. Lett. 64, 127 (1990). 

E. Brezin and V. A. Kazakov, Phys. Lett. B 236, 144 (1990). 

D. J. Gross and N. Miljkovic, Phys. Lett. B 238, 217 (1990). 

E. Brezin, V. A. Kazakov and A. B. Zamolodchikov, Nucl. Phys. B 338, 673 (1990). 
P. H. Ginsparg and J. Zinn-Justin, Phys. Lett. B 240, 333 (1990). 

H. Dorn and H. J. Otto, |arXiv:hep-th/9501019| 

A. B. Zamolodchikov and A. B. Zamolodchikov, Nucl. Phys. B 477, 577 (1996) 
|arXiv:hep-th/9506136| . 

J. Teschner, Phys. Lett. B 363, 65 (1995) |arXiv:hep-th/9507109| . 
V. Fateev, A. B. Zamolodchikov and A. B. Zamolodchikov, arXiv:hep-th/0001012 



J. Teschner, arXiv:hep-th /0009138 

J. Teschner, Class. Quant. Grav. 18, R153 (2001) |arXiv:hep-th/0104158| . 



A. B. Zamolodchikov and A. B. Zamolodchikov, arXiv:hep-th/0101152 , 

B. Ponsot and J. Teschner, Nucl. Phys. B 622, 309 (2002) |arXiv:hep-th/0110244| . 

T. Fukuda and K. Hosomichi, Nucl. Phys. B 635, 215 (2002) |arXiv:hep-th"7"0 202032 ; 

C. Ahn, C. Rim and M. Stanishkov, Nucl. Phys. B 636, 497 (2002) 
|arXiv:hep-th /0202043] . 

J. McGreevy and H. Verlinde, JHEP 0312, 054 (2003) |arXiv:hep-th/0304224| . 

I. R. Klebanov, J. Maldacena and N. Seiberg, JHEP 0307, 045 (2003) 
|arXiv:hep-th/0305159| . 

T. Takayanagi and N. Toumbas, JHEP 0307, 064 (2003) |arXiv:hep-th/0307083| . 

M. R. Douglas, I. R. Klebanov, D. Kutasov, J. Maldacena, E. Martinec and N. Seiberg, 
arXiv:hep-th/0307195 



[19] I. R. Klebanov, J. Maldacena and N. Seiberg, arXiv:hep-th/0309168 



31 



[20] E. J. Martinec, |arXiv:hep-th /0305148 , 

[21] J. McGreevy, J. Teschner and H. Verlinde, JHEP 0401, 039 (2004) 
|arXiv:hep-th/0305194| . 

[22] S. Y. Alexandrov, V. A. Kazakov and D. Kutasov, JHEP 0309, 057 (2003) 
|arXiv:hep-th/0306T77| . 



[23] A. Sen, Mod. Phys. Lett. A 19, 841 (2004) arXiv:hep-th/0308068 ; JHEP 0405, 076 
(2004) |arXiv:hep-th/0402157| . |arXiv:hep-th/0408064| 

[24] J. McGreevy, S. Murthy and H. Verlinde, JHEP 0404, 015 (2004) 
|arXiv:hep-th/0308105l . 

[25] J. L. Karczmarek and A. Strominger, JHEP 0404, 055 (2004) |arXiv:hep-th7 0309138 . 

[26] O. DeWolfe, R. Roiban, M. Spradlin, A. Volovich and J. Walcher, JHEP 0311, 012 
(2003) arXiv:hep-th/0309148|. 



[27] S. Alexandrov, arXiv:hep-th/0310135 



[28] J. Gomis and A. Kapustin, JHEP 0406, 002 (2004) |arXiv:hep-th/0310195| . 

[29] C. V. Johnson, JHEP 0403, 041 (2004) |arXiv:hep-th/0311129| . 

[30] D. J. Gross and J. Walcher, JHEP 0406, 043 (2004) |arXiv:hep-th/0312021| . 

[31] M. Aganagic, R. Dijkgraaf, A. Klemm, M. Marino and C. Vafa, |arXiv:hep-th/ 0312085 . 

[32] N. Seiberg and D. Shih, JHEP 0402, 021 (2004) |arXiv:hep-th/0312170| . 



[33] J. de Boer, A. Sinkovics, E. Verlinde and J. T. Yee, JHEP 0403, 023 (2004) 
|arXiv:hep-th/0312T35| . 



[34] D. Gaiotto and L. Rastelli, |arXiv:hep -th/0312196 . 



[35] I. K. Kostov, Nucl. Phys. B 689, 3 (2004) |arXiv:hep-th/0312301| . 
[36] D. Berenstein, JHEP 0407, 018 (2004) |arXiv:hep-th/0403110| . 
[37] V. A. Kazakov and I. K. Kostov, |arXiv:hep- th/0403152 , 



[38] A. Kapustin and A. Murugan, arXiv:hep-th/0404238 

[39] M. Hanada, M. Hayakawa, N. Ishibashi, H. Kawai, T. Kuroki, Y. Matsuo and T. Tada, 
Prog. Theor. Phys. 112, 131 (2004) |arXiv:hep-th/0405076| . 

32 



[40] D. Kutasov, K. Okuyama, J. w. Park, N. Seiberg and D. Shih, JHEP 0408, 026 (2004) 
|arXiv:hep-th/0406030| . 



[41 
[42 
[43; 
[44 
[45 
[46; 
[47; 
[48; 
[49 
[50 

[51 

[52; 

[53; 
[54 
[55 
[56; 
[57; 

[58; 



D Ghoshal, S. Mukhi and S. Murthy, |aIXiv7he p-th/0406106 



J. Ambjorn, S. Arianos, J. A. Gesser and S. Kawamoto, arXiv:hep-th/0406108 
E. Martinec and K. Okuyama, |arXiv:hep-tn7 0407136 . 



J. Maldacena, G. W. Moore, N. Seiberg and D. Shih, |arXiv:hep-t h/0408039 . 

C. V. Johnson, |ar"Xiv:hep-th/0408049[ 

J. E. Carlisle and C. V. Johnson, arXiv:hep-th/0408159 



N. Itzhaki and J. McGreevy, |arXiv:hep-th /0408180 

S. Giusto and C. Imbimbo, |arXiv:hep-th7 0408216 

D. J. Gross and E. Witten, Phys. Rev. D 21, 446 (1980). 

V. Periwal and D. Shevitz, Phys. Rev. Lett. 64, 1326 (1990); V. Periwal and 
D. Shevitz, Nucl. Phys. B 344, 731 (1990). 

C. Crnkovic, M. R. Douglas and G. W. Moore, Nucl. Phys. B 360, 507 (1991); 

C. Crnkovic, M. R. Douglas and G. W. Moore, Int. J. Mod. Phys. A 7, 7693 (1992) 
|arXiv:hep-th/9108014| . 

R. C. Brower, N. Deo, S. Jain and C. I. Tan, Nucl. Phys. B 405, 166 (1993) 
|arXiv:hep-th/9212127| . 

T. J. Hollowood, L. Miramontes, A. Pasquinucci and C. Nappi, Nucl. Phys. B 373, 247 
(1992) |arXiv:hep-th/9109046| . 

S. Dalley, C. V. Johnson and T. R. Morris, Nucl. Phys. B 368, 625 (1992); S. Dalley, 
C. V. Johnson, T. R. Morris and A. Watterstam, Mod. Phys. Lett. A 7, 2753 (1992). 

R. Lafrance and R. C. Myers, Mod. Phys. Lett. A 9, 101 (1994) 
|arXiv:hep-th/9308113| . 

I. K. Rostov, Mod. Phys. Lett. A 4, 217 (1989). 

S. Dalley, C. V. Johnson and T. R. Morris, Nucl. Phys. B 368, 625 (1992). Nucl. Phys. 
B 368, 655 (1992). 



33 



[59] V. A. Kazakov, A. A. Migdal and I. K. Rostov, Phys. Lett. B 157, 295 (1985); 

D. V. Boulatov, V. A. Kazakov, I. K. Kostov and A. A. Migdal, Nucl. Phys. B 275, 
641 (1986). 



[60 
[61 
[62 
[63 
[64 
[65 
[66 

[67; 
[68; 

[69 

[70; 

[71 
[72 
[73; 
[74 

[75; 

[76; 

[77; 
[78; 



F. David, Phys. Lett. B 159, 303 (1985). 

I. K. Kostov and M. L. Mehta, Phys. Lett. B 189, 118 (1987). 
I. R. Klebanov and R. B. Wilkinson, Phys. Lett. B 251, 379 (1990); 
I. R. Klebanov and R. B. Wilkinson, Nucl. Phys. B 354, 475 (1991). 
J. D. Edwards and I. R. Klebanov, Mod. Phys. Lett. A 6, 2901 (1991). 
J. C. Plefka, arXiv:hep-th/9601041 



H. G. Kausch, |arXiv:hep-th/9510149| 

H. G. Kausch, Nucl. Phys. B 583, 513 (2000) |arXiv:hep-th/0003029| . 

M. A. I. Flohr, Int. J. Mod. Phys. A 12, 1943 (1997) |arXiv:hep-th/9605151| . 

M. R. Gaberdiel and H. G. Kausch, Phys. Lett. B 386, 131 (1996) 
|arXiv:hep-t"*h /9606050|. 

M. R. Gaberdiel and H. G. Kausch, Nucl. Phys. B 538, 631 (1999) 
|arXiv:hep-th/9807091| . 

I. I. Kogan and J. F. Wheater, Phys. Lett. B 486, 353 (2000) |arXiv:hep-th/0003184 



S. Kawai and J. F. Wheater, Phys. Lett. B 508, 203 (2001) arXiv:hep-th/0103197 . 



A. Bredthauer and M. Flohr, Nucl. Phys. B 639, 450 (2002) arXiv:hep-th/0204154 . 

P. Di Francesco, P. H. Ginsparg and J. Zinn-Justin, Phys. Rept. 254, 1 (1995) 
|arXiv:hep-th/9306153| . 

G. Bonnet, F. David and B. Eynard, J. Phys. A 33, 6739 (2000) 
|arXiv:cond-mat/000332"4| . 

R. Dijkgraaf, S. Gukov, V. A. Kazakov and C. Vafa, Phys. Rev. D 68, 045007 (2003) 
|arXiv:hep-th/0210238| . 

B. Duplantier and I. K. Kostov, Nucl. Phys. B 340, 491 (1990). 
I. K. Kostov, Phys. Lett. B 266, 42 (1991). 

34 



[79] I. K. Rostov and M. Staudacher, Nucl. Phys. B 384, 459 (1992) 
|arXiv:hep-th/9203030) . 



[so; 

[81 
[82 

[83; 
[84 

[85 

[86; 
[87; 

[88 

[89; 

[90 
[91 
[92 

[93 

[94 
[95; 

[96 
[97 



B. Eynard and J. Zinn-Justin, Nucl. Phys. B 386, 558 (1992) |arXiv:hep-th7 9204082 . 

B. Eynard and C. Kristjansen, Nucl. Phys. B 455, 577 (1995) |arXiv:hep-th7 9506193 . 

V. A. Kazakov and I. K. Rostov, Nucl. Phys. B 386, 520 (1992) 
|arXiv:hep-th/9205059| . 

I. R. Rostov, Nucl. Phys. B 658, 397 (2003) |arXiv:hep-th/0212194| . 

I. R. Rostov, B. Ponsot and D. Serban, Nucl. Phys. B 683, 309 (2004) 
|arXiv:hep-th/0307189| . 

J. M. F. Labastida, M. Pernici and E. Witten, Nucl. Phys. B 310, 611 (1988). 
E. Witten, Surveys DifT. Geom. 1, 243 (1991). 
J. Distler, Nucl. Phys. B 342, 523 (1990). 

R. Dijkgraaf, H. Verlinde and E. Verlinde, Nucl. Phys. B 348, 435 (1991). 
E. Witten, Nucl. Phys. B 268, 253 (1986). 

A. LeClair, M. E. Peskin and C. R. Preitschopf, Nucl. Phys. B 317, 411 (1989). 
M. R. Douglas and G. W. Moore, |arX"iv:hep-th/9603167 



V. G. Rnizhnik, A. M. Polyakov and A. B. Zamolodchikov, Mod. Phys. Lett. A 3, 819 
(1988). 

F. David, Mod. Phys. Lett. A 3, 1651 (1988); J. Distler and H. Rawai, Nucl. Phys. B 
321, 509 (1989). 

E. J. Martinec, G. W. Moore and N. Seiberg, Phys. Lett. B 263, 190 (1991). 

S. Govindarajan, T. Jayaraman and V. John, Nucl. Phys. B 402, 118 (1993) 
|arXiv:hep-th/9207109| ; Int. J. Mod. Phys. A 10, 477 (1995) |arXiv:hep-th/9309040| . 

J. L. Cardy, J. Phys. A 25, L201 (1992) |arXiv:hep-th/9111026| . 

J. S. Caux, I. Rogan, A. Lewis and A. M. Tsvelik, Nucl. Phys. B 489, 469 (1997) 
|arXiv:hep-th /9606138 . I. I. Rogan and A. Nichols, Int. J. Mod. Phys. A 17, 2615 
(2002) |arXiv:hep-th/0107160| . 



35 



[98] B. H. Lian and G. J. Zuckerman, Phys. Lett. B 254, 417 (1991). 

[99] M. Bershadsky, W. Lerche, D. Nemeschansky and N. P. Warner, Nucl. Phys. B 401, 
304 (1993) arXiv:hep-th/9211040| . 

[100] A. Losev, Theor. Math. Phys. 95, 595 (1993) [Teor. Mat. Fiz. 95, 307 (1993)] 
|arXiv:hep-th/92 1 1090) . 



[101] W. Lerche, arXiv:hep-th/9401121 



[102] M. Bershadsky and I. R. Klebanov, Phys. Rev. Lett. 65, 3088 (1990); Nucl. Phys. B 
360, 559 (1991). 

[103] D. V. Boulatov and V. S. Kazakov, Phys. Lett. B186 (1987), 379. 



36 



